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Mixtures of polarised fermions of two different masses can form weakly-bound clusters, such 
as dimers and trimers, that are universally described by the scattering length between the heavy 
and light fermions. We use the resonating group method to investigate the low-energy scattering 
processes involving dimers or trimers. The method reproduces approximately the known particle- 
dimer and dimer-dimer scattering lengths. We use it to estimate the trimer-trimer scattering length, 
which is presently unknown, and find it to be positive. 


I. INTRODUCTION 

In the last decade, the use of controlled Feshbach 
resonances in ultra-cold atom experiments have en¬ 
abled the study of low-energy quantum systems of par¬ 
ticles interacting with large scattering lengths. Close 
to a Feshbach resonance, the interparticle scattering 
length is much larger than the range of interparti¬ 
cle forces. As a result, the low-energy properties of 
these systems are universal, in the sense that they 
depend only upon a few parameters, such as the scat¬ 
tering length [T], and the three-body parameter Eiia 
in systems exhibiting the Efimov effect @H6] . More¬ 
over, close to Feshbach resonances, atoms can be as¬ 
sociated into clusters of universal character: diatomic 
molecules called Feshbach molecules that are a reali¬ 
sation of universal dimers mm, triatomic molecules 
that are a realisation of Efimov states mm- The¬ 
ory predicts the existence of a variety of other univer¬ 
sal clusters of larger number of particles jTBM20] that 
are expected to be observed experimentally in the fu¬ 
ture EJ. 

The few-body properties, in particular the scatter¬ 
ing properties of clusters, can play a crucial role in the 
identification and stability of the many-body ground 
states of these systems. For instance, the stability of a 
gas of universal dimers made of fermions was observed 
0HMU and explained theoretically f22] E3| by exact 
four-body calculations for two scattering dimers. 

Although it is sometimes feasible to calculate ex¬ 
actly the wave function of an IV-body cluster [JJ [23], 
the exact computation of the scattering properties of 
two clusters is generally out of reach for N > 3. In the 
context of nuclear and sub-nuclear physics, a broad ar¬ 
ray of approximation schemes have been successfully 
developed to address similar problems. One of the 
leading techniques is the so-called Resonating Group 
Method (RGM), introduced by Wheeler [25], to study 
light nuclei, such as 16 O and 8 Be, modelled as clusters 
of a particles. Since then, it has been employed in a 
variety problems including the scattering of light nu¬ 
clei, the stability of light nuclei to external nucleon 
scattering and nuclear particles EH E3- More re¬ 
cently, it has been used |28U30] to study low-energy 
scattering, and bound states, of baryon-baryon and 


other multi-quark cluster configurations. 

In the single-channel approximation, the RGM con¬ 
structs the low-energy scattering wave function of two 
or more scattering clusters from the wave functions of 
the individual clusters, while preserving the full anti- 
symmetrization of wave functions. This gives an effec¬ 
tive potential between the clusters that can be used to 
treat scattering as well as bound states. It is especially 
accurate in situations in which single clusters are not 
strongly altered by the scattering process. Here, we 
propose to apply this method to the low-energy scat¬ 
tering of universal fermionic clusters that are relevant 
to ultra-cold atoms close to Feshbach resonances. 

The paper is organised as follows. In section [IT] we 
review the essence of the RGM. In section |III| we ap¬ 
ply it to universal clusters whose scattering properties 
are known. In section |IV| we apply the RGM to the 
yet unknown scattering of universal trimers. 

II. THE RESONATING GROUP METHOD 
A. General formalism 

Let us consider the scattering between a cluster A 
of n particles and a cluster B of N — n particles. It is 
assumed that the wave functions <^ A (1, 2,..., n) and 
<t>n{n + 1, n + 2,..., N) of these clusters are known. 
In the single-channel RGM, the 7V-body wavefunction 

describing the scattering process is constructed as 
the antisymmetrised product of the cluster wave func¬ 
tions and a wave function V’(-R) for the relative motion 
between the two clusters: 

A/ =S [0 A (1, 2,..., n)(/) B (n + 1, n + 2,..., N)if>(R)] . 

(!) 

Here, S denotes the symmetrisation (or antisym- 
metrisation) operator that symmetrises (or antisym¬ 
metrises) the wave function under the exchange of 
identical particles. Symmetrisation is performed for 
bosonic particles, whereas antisymmetrisation is per¬ 
formed for fermionic particles. The vector R describes 
the relative position between the centres of mass of the 
two clusters. The idea behind this approximation is 
that the structure of the two clusters is not much al¬ 
tered during the collision, and the two clusters mix 
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only through the exchange of identical particles. 

The purpose of the RGM is to determine the wave 
function ip{R) for the relative motion of the clusters. 
This is done by applying the variation principle to the 
average quantity, 


(2) 

where H is the total hamiltonian and £ is the total 
energy of the system. Requiring ip to extremise the 
above quantity implies that for an infinitesimal varia¬ 
tion Sip around ip we have 


where E = £ — E A — E B is the scattering energy be¬ 
tween the two clusters. The symmetrisation operator 
S can be written as S = 1 + 5', i.e. the action of S 
gives one term leaving the wave function unchanged, 
and other terms where particles are exchanged. Thus, 
Eq. § can be written as 


(1 - K) • (' T r - E)ip + V B iP + V BX1 -iP = 0 


( 10 ) 


where we have introduced a local potential 14 called 
the direct potential, 


(S[(p A (p B 6'ip\\H - £\S[<p A (p B ip]) +c.c. = 0. 


V d (R) = {(p A (p B \V AB \(p A (p B ), (11) 


The variations Sip and its complex conjugate Sip* can 
be formally taken as independent variations, resulting 
in the following Euler-Lagrange equation of motion, 

{S[cp A (p B )\H - £\S[cp A (p B ip]) = 0, (3) 

which can be simplified as 

(cp A (p B \H - £\S[(p A cp B 'ip]) = 0. (4) 

since the total hamiltonian H is invariant under the 
exchange of identical particles. 

The hamiltonian H consists of kinetic operators tj 
for each particle and pairwise interaction terms Vij for 
each pair of particles, 

N N 

+ ZX’ - *C- (5) 

i —1 i<j 

We have subtracted the kinetic operator t c for the 
centre of mass, since it can be eliminated from the 
problem. The hamiltonian can be rewritten as 

H = H a + H b +T r + V ab , (6) 

where H A and H B denote the internal hamiltonian of 
each cluster A and B, Tr denotes the kinetic operator 
for the relative motion of the two clusters, and V AB is 
the sum of interactions between the two clusters. The 
wave functions cp A and <p B are eigenstates of H A and 
H b with eigenvalues E A and E B , i.e. 

H A (p A = E A cp A and H B (p B = E B (p B . (7) 

There are two ways this can be used to simplify the 
equation of motion Eq. Q. Either one applies the 
hamiltonian to the wave functions <f> A and </> B on the 
right-hand side, or to the wave functions cp A and cp B 
on the left-hand side. We refer to these two equivalent 
procedures as the RGM1 and RGM2. Although they 
result in formally different equations, their solutions 
are the same. 

In the RGM1, one writes 

(4> A (p B \S[(H - £)(p A (p B ip]) = 0 (8) 

and using Eqs. © and 0 

(cp A cp B \S[{T R + H AB - E)<f>M]) = 0, (9) 


a non-local potential V EX1 called the exchange poten¬ 
tial, 


V EX1 -ip = J d 3 R'V EX1 (R,R')iP(R') (12) 

= (cpA<pB\S'[V AB (p A <p B ip]), 

and a non-local operator K called the exchange kernel, 

K ■ ip = J d 3 R'K{R , R')ip(R') 

= -{cp A (p B \S'[(p A (p B ip]}. 


(13) 


In the RGM2, one applies the hamiltonian Eq. ^ 
to the wave functions <p A (p B on the left-hand side of 
Eq. using Eq. (J7|. This gives 

<<Mb|C z h + ^ab - E)S[(p A <p B ip]) = 0, (14) 

which can be written as 

{Tr - E)( 1 - K) • iP + V B iP + I4 X2 • iP = 0 (15) 

where the exchange potential Ve X 2 is defined by 

14x2 •V’ = J d 3 R'V EX2 {R, R^ipiR 1 ) 

= (<p A <p B \V AB S'[cP A cP B iP]). (16) 

Hence, the RGM consists in calculating the poten¬ 
tials 14, 14 x and kernel K , and solving the equa¬ 
tion for the relative motion between the two clus¬ 
ters, either Eq. (10) or (151. This is of course a great 
simplification over solving the full iV-body equation. 
Nonetheless, the determination of 14, 14x and K in¬ 
volve 3(n — 1) + 3(IV — n — 1) = 3{N — 2 (-dimensional 
integrals whose computation may be costly for large 
N. 


B. RGM with contact interactions 

In the following, we apply the RGM to the scatter¬ 
ing of universal clusters. Their universal character is 
described by the zero-range theory, which corresponds 
to the limit of the range of interaction being much 
smaller than the s-wave scattering length a. In this 
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limit, the interaction potential V\ 3 between two parti¬ 
cles appearing in Eq. (|5| and included in the term V AB 
in Eq. ([6]) can be approximated by a contact potential, 

Vij(r) = g6 3 {r)-^-r- (17) 

with the coupling constant 


9 = 


4:nh 2 a 

2/i 


(18) 


Here, fi is the reduced mass of the two interacting par¬ 
ticles, and d/drr■ is an operator regularising the 1/r 
divergence of the wave function when particles come 
into contact (r = 0). This potential binds two parti¬ 
cles only for a > 0, and we restrict our consideration 
to this case throughout this paper. 

The presence of the three-dimensional Dirac delta 
function in the potential Eq. reduces by three the 
dimensionality of the integrals. The dimensionality of 
the integrals Eqs. (Ill, ( fl2j ), and ( jl6| , for V D , V EX1 
and 14x2, is thus reduced to 3(TV — 3)7 


C. Partial-wave expansion 

To proceed further, one can perform a partial-wave 
expansion in spherical harmonics Y)m i n the RGM1 
and RGM2 equations. The relative wave function is 
expanded as 

'ip(R) = E tfl>tm(R)Y tm (R), (19) 

im 


where R denotes the orientation of R. Then, the 
RGM1 equation, Eq. (101, becomes the following set 
of coupled equations: 

(Tk - E) 4>tm(R) + E < m/m ' ( R ) (20) 

t'm' 

pOO 

- V / dR' K trn ’ t ' rn ' (R, R') (T& - E)ip Vm ,{R!) 

I'm ' J ° 

/»oo 

+ E / dR' V^/ m '(R,R')^m'(R') = 0 

Dlrr,/ ^0 


and the RGM2 equation, Eq. (15), becomes the set of 
coupled equations 


(Tk ~ E) t em (R ) + E Vi m ' l ' m \R^ Vm ,{R) (21) 

£'m' 

roo 

-E / dR/ (Tk-E)K tm * m '{R,R!)il> l . rn ,(R!) 

I'm' J ° 

pOO 

+ E / dRl V^/ m \R,R')^ m '(R') = 0 

I'm' 


where /j,n is the reduced mass of the two clusters, and 
V‘ m ’ t ' m \R) = jd 2 R Y e * m (R)V D (R)Y rm ,(R), (23) 


K im,e'm' = RR' Jd 2 Rd 2 R' (24) 

Y; m (R)K{R,Rf)Y em ,{Rf), 


V„ 


£m£ m 


{R, R’) = RR' d 2 Rd 2 R’ 


(25) 


Y; m (R)V EX (R,R')Y eim ^R'). 


The dimensionality of integration in Eqs. (231, 


and (251 is, respectively, 3(TV — 2) — 1, 3 (N — 2) + 1, 
and 3 JN - 2) - 2. 


D. Local approximation 


It turns out, as we shall see in the cases treated 
below, that the contribution from the non-local kernel 
K is often small and may be neglected. Moreover, 
in some cases, the exchange potentials m and 

V,fx 2 d m are nearly local and may be approximated 
by the local potentials 


R E x m £ir'(R)= / dR'V^'m' {R R/) (26) 

JO 


Neglecting K and using the local form Eq. (261 


of the exchange potentials constitute the local RGM 
approximation. In this approximation, RGM1 and 
RGM2 equations have the form of conventional 
Schrodinger equations: 


(Tk ~ E) il> tm (R) + E Vf™’*' 171 '(R)i/>t>m'(R) = 0 

i'm' 

(27) 


(Tk - E) 1> tm {R) + E Vimi'm'(R) = 0 


, (28) 

i t rZm.i m t rtm.t m , Trim.cm' i T rim.t m 

where v\ — vE H - Vexi local 5 v 2 

| t r£m,£'m' 

D I 'EX 2 local • 


Unlike the RGM1 and RGM2 equations, Eqs. (201 
and (211, the local RGM1 and RGM2 equations, 
Eq. (271 and (281, are not equivalent. Nevertheless, 
they often yield similar results as we shall see in the 
following sections. 


E. Scattering length and scattering volume 


with the kinetic energy operator 

h 2 ( d 2 £[t- 


n = 


2^ n 


d 2 

dR 2 


1) 


R 2 


( 22 ) 


After solving the RGM equations in partial waves, 
Eq. (20) o r fl2l] ), or their local-potential approxima¬ 
tion, Eq. (|27[) or (28), one obtains the partial-wave 
components ijjtm(R) of the relative wave function if). 
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Figure 1: Schematic representation of the two permuta¬ 
tions of identical fermions between a dimer and a fermion. 


For zero scattering energy (E = 0), one can extract 
the partial-wave scattering lengths from these compo¬ 
nents. 

From the s-wave component V’oo (R) ex i? + 0(l) for 
R —> oo, one obtains the s-wave scattering length, 


a = lim R — 
R—>oo 


tl>oo(R) 


(29) 


are two possibilities, as shown in Fig. |T] no permuta¬ 
tion and the exchange of two fermions of mass M. It 
follows that the direct and exchange potentials of the 
RGM equations, Eqs. (101, and (151, are given by the 
following expressions: 


V D (R) = g 


1 


vi —R 


(32) 


F EX1 • ^R) = gg>*(0)<p(R)1>( - ^ R ) (33) 

Wfi) = w(0) 

(34) 

where k = M/m is the mass ratio and 

£(0) = - I™ S~r ■ ip{f) = J— > ( 35 ) 

r->o dr \Z2naa 

The exchange kernel is given by 


and from the p -wave component Mm(R) oc R 2 
0(1/R), one obtains the p -wave scattering volume, 


v = lim 


1 M R 3 RMlm(R) - 2 Mm(R) 

3 RMim(R) + i’lm(R) ' 


(30) 


These formulas follow from the standard definition of 
the scattering phase shifts ED- 


III. SCATTERING OF UNIVERSAL DIMERS 
A. Universal dimers 


k ■ MR) = 


(« + i ) 6 

(1 + 2 k ) 3 


d 3 R'ip*(r i)(^(r 2 )V’(-R') 


with 


(36) 


r i = 


r 2 = 


(^ + 1) S 
1 -f- 2 K 
k{k + 1 ) 

1 + 2 K 


R' 

R! 


k(k + 1 ) 
1 + 2k 

(* + l) 2 

1 H- 2ac 


R 

R. 


The kinetic operator in Eqs. (101, (151 is given by 

h 2 

Tr = - 


2 ( 


M+m 


-xpi v «- 

M > 


We consider universal dimers made of a polarised 
fermion of mass M and a polarised fermion of mass 
to. These dimers are two-body s-wave weakly-bound 
states. The normalised wave function <p( r ) f° r the 
relative motion of the two particles inside the dimer 
is given by 


The RGM1 and RGM2 equations can be solved by 
performing the partial-wave expansion of section IIC 


Here, the potentials Eqs. (32p4l do not couple partial 


waves: 


t rim,£'m' 


CK $ra,rn/ 


(37) 


Mr) = 



(31) 


B. Scattering of a dimer and a particle 

First, we consider the scattering of a universal 
dimer with a fermionic particle of mass M. To ap¬ 
ply the RGM to this case, we set </> A = ip given by 
Eq. @ (f> B = 1, and the interaction potential given 
by Eq. ©, assuming that there is no interaction be¬ 
tween identical fermions. The antisymmetrisation op¬ 
erator in the calculation of the exchange potentials 
and kernel is obtained by considering all possible per¬ 
mutations of identical fermions. In this case, there 


V 'Um.i'm' e r 

ex ^ Vm,n 


(38) 


and for a given partial wave (t, to), we obtain from 
Eq. p5l, 


V^/ m (R,R f ) = (-l) e g—M(0) 


(39) 


xm r )s(r , - 1 ^- 1 r), 


V^/ m (R,R') = (-1 Yg 


K - I - 1 


£(0) (40) 


xM(^R)s(r' - ^r). 


The factor (—l) f in these expressions comes from 
the minus sign in the argument of 0 hr Eqs. (33p4 1. 
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Figure 2: Fermion-dimer s-wave scattering length a/d in 
units of a as a function of the mass ratio. 



Figure 3: Fermion-dimer p-wave scattering volume Vfd in 
units of a 3 as a function of the mass ratio. 


Due to this factor, the exchange potential is repulsive 
for even partial waves, and it is attractive for odd par¬ 
tial waves. Moreover, Eqs. ( 39]|40 l show that the ex¬ 
change potentials have an increasingly local character 
as the mass ratio k increases. Their local approxima¬ 


tion, given by Eq. (26), leads to 


= (-i Yg 


K + 1 


(OMR), 


(41) 


= (-1 Yg 


K + 1 




We solve the resulting RGM and local RGM equa¬ 
tions numerically by discretising the coordinate R. 


1. s-wave scattering 

We first consider fermion-dimer scattering in the 
s wave, for which the effective potential is purely re¬ 
pulsive. The fermion-dimer s-wave scattering length 
ajd is therefore always positive. It is shown in Fig. [2] 
as a function of the mass ratio k. For the equal mass 
case (M = to), we obtain 


afd = 1.19 a, 


2. p-wave scattering 


In the p-wave channel, the fermion and dimer at¬ 
tract each other. This is due to the Efimov at¬ 
traction 0 0 that results from the effective inter¬ 
action between the two heavy fermions mediated by 
the light fermion. Although the Efimov attraction 
wins over the centrifugal barrier only for mass ra¬ 
tios M/m > k c ~ 13.6069657 [6J [2H 31], resulting 
in an infinite discrete-scale-invariant tower of three- 
body bound states, it also makes the system attrac¬ 
tive for lower mass ratios, resulting in an overall neg¬ 
ative p-wave scattering length. As the mass ratio in¬ 
creases, the Efimov attraction strengthens, and two 
universal three-body bound states appear at mass ra¬ 
tios K\ — 8.17260 and K 2 = 12.91743 [21]. At these 
mass ratios, fermion-dimer p -wave scattering is reso¬ 
nant and the p-wave scattering volume vjd diverges, 
as shown in Fig. [3] 

In the RGM, the effective potential between the 
fermion and the dimer scattering in the p wav e is als o 
attractive, due to the factor (—1)^ of Eqs. ( 39p0 1. 
The scattering volume calculated in the RGM is thus 
negative and very close to the exact one up to the mass 
ratio M/m ss 6. For the equal mass case (M = to), 
the RGM gives 


which is consistent with the exact result ss 
1.17907a [32 ] G2J. All RGM results are within 2% 
of the exact results, indicating that there is little exci¬ 
tation during the collision of a dimer and fermion, the 
dimer remaining bound during the collision. Nonethe¬ 
less, the exchange of particles is crucial. The dotted 
curve in Fig. [2] shows that including only the direct 
potential (neglecting the exchange kernel and poten¬ 
tial) yields a much smaller scattering length. On the 
other hand, the exchange kernel K brings a significant 
difference only for mass ratios smaller than one, and 
may be neglected otherwise, as shown by the dashed 
red and blue curves in Fig. [2] As to the local approx¬ 
imation, it leads to results which are close to those 
of the RGM for sufficiently large mass ratios, as seen 
from the red and blue curves in Fig. [2j 


Vfd = —0.98a, 

which is consistent with the exact result ss 
—0.96a [31]. Beyond the mass ratio ~ 6, the RGM 
results deviate strongly from the exact results. This 
is explained by the fact that the resonance and the 
three-body bound state at M/m = Ki imply three- 
body correlations that are not fully captured by the 
RGM. Nevertheless, the RGM exhibits a similar reso¬ 
nance, but at a shifted mass ratio k^ rgm) « 9.5. This 
shows that the Efimov attraction, physically due to 
the exchange of light fermion between the two heavy 
fermions, is partially captured by the mere antisym- 
metrisation of the wave function in the RGM, as sug¬ 
gested by Fig. |T] 

The local RGM equations reproduce approximately 
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Figure 4: Schematic representation of the two permuta¬ 
tions of identical fermions between two dimers. 


the RGM results for M/m < 6, as shown by the blue 
and red curves in Fig. [3] For larger mass ratios, the 
difference between the local RGM and full RGM re¬ 
sults is substantial and it is mainly due to the absence 
of the exchange kernel K in the local RGM equations, 
as shown by the dashed red and blue curves in Fig. [3] 


C. Scattering of two dimers 



0 20 40 60 80 100 


Mass ratio k 

Figure 5: Dimer-dimer scattering length add in units of a 
as a function of the mass ratio. 


The kinetic operator in Eqs. (101, (151, (27) and 


(28) is given by 


Now, we consider the scattering of two universal 
dimers. We thus apply the RGM equations for the two 
cluster wave functions <p A = <f> B = ip given by Eq. (311 


and the interaction potential given by Eq. (171, assum¬ 


ing again that there is no interaction between identi¬ 
cal fermions. The antisymmetrisation operator in the 
calculation of the exchange potentials and kernel is 
obtained by considering all possible permutations of 
identical fermions. In this case, there are two possi¬ 
bilities, as shown in Fig. [4j no permutation and the 
exchange of two fermions of mass M (which is equiva¬ 
lent to exchanging the two fermions of mass m). After 
some straightforward calculations, the direct and ex¬ 
change potentials, as well as the exchange kernel of 
the RGM equations, Eq. (101 and (15), are given by 
the following expressions: 

V D (R) =2 g(n + l) 3 J d 3 r\ip(r)ip((K + 1)R + kt-)| 2 , 


K(R , R!) = — A/ d 3 rip*(r+R 3 ,)ip*(r)(p(r+Ri)ip(r+R 2 ), 

V EX1 (R,R') = 2g\ip* (R 3 )(p* (0)ip(Ri)ip(R 2 ), 

V EX2 (R,R') = 2g\ip*(R 2 )<p*(R 1 )v(0)cp(R 4 ). 

In these expressions, we have set n = M/m, A = 
(k + l) 6 /(2 k) 3 , and 

R i = ^(R' + R), 

R 2 = ^l(R'-R), 

R 3 = {k ^ I ((k + 1)R'-(k-1)R), 

Ri = ((k - 1)R' - (k + l)R ). 


T r 


h 2 

M + m 


V 2 


R- 


(43) 


We solve the RGMl and RGM2 equations, as well as 
their local approximation, by performing the partial- 
wave expansion of section |II C| and discretising the co¬ 
ordinate R. The potentials are repulsive for all partial 
waves. The resulting dimer-dimer s-wave scattering 
length is shown in Fig. [5] as a function of the mass 
ratio n. For the equal mass case (M = m), we obtain 


add = 0.752 a, 


which is close, although significantly different, from 
the exact result ss 0.6a |22l 23] [3HH3T] . This means 
that compared to dimer-particle scattering, there is a 
bit more excitation during the collision of two dimers, 
although it remains small. On the other hand, ex¬ 
change is less important than in the case of fermion- 
dimer scattering, as the major contribution to the 
scattering length comes from the direct potential, as 
seen from the dotted curve in Fig. [5] The exchange 
potentials have an increasingly local character as the 
mass ratio increases. We have calculated the scatter¬ 
ing length with the RGM up to mass ratio 20. Beyond 
this mass ratio, the local character of the potential 
makes it difficult to solve the problem as a non-local 
one, since a high degree of discretisation is needed. 
The local RGM equations, on the other hand, are eas¬ 
ier to solve. They give results which are very close to 
the RGM, as can be seen from the blue and red curves 
of Fig. [5] and can easily be extended to larger mass 
ratios. 



























7 


© O 

R 

r 


Figure 6: Jacobi coordinates r and R decribing a trimer 
made of two heavy fermions and a light fermion. The vec¬ 
tor r = r 2 — ri is the relative position between particle 
2 and 1, and the vector R = r 3 is the rela¬ 

tive position between particle 3 and the centre of mass of 
particles 1 and 2. 


IV. SCATTERING OF UNIVERSAL 
TRIMERS 


A. Universal trimers 

We now consider universal trimers made of two po¬ 
larised fermions of mass M and a polarised fermion 
of mass to. Such trimers exist for a mass ratio 
M/m > Ki ~ 8.17260. They rotate with one quan¬ 
tum unit of angular momentum, and can therefore be 
in three possible internal quantum states of rotation, 
labelled by the quantum number to € {—1,0,1}. For 
a mass ratio M/m > k c ~ 13.6069657 [BJ Ml M], the 
trimers are Efimov states mu], characterised by the 
scattering length a between the two different kinds of 
fermions, and a three-body parameter. For a mass 
ratio M/m < n c , the trimers are Kartavtsev-Malykh 
states |M], characterised only by the scattering length 
a. We restrict our consideration to these states, and 
therefore to the range k± < M/m < k c where a 
ground-state trimer exists. 

The trimer wave function is expressed as a func¬ 
tion of Jacobi vectors r and R shown in Fig. [u] To 
a good accuracy, the trimer wave function is well ap¬ 
proximated by the adiabatic hyperspherical form [2?i 


*(r,R) = 


W) 


H 5 / 2 [ 

where the component 


i/j Fad {TZ 1 a,y) - i/ 

Fad (11, a, y) 


(44) 


i/j Ftld (TZ,a,y) 


cm 

sin 2 a 


Vw( n,a)YT{y) 


(45) 


incorporates the angular momentum of the trimer 
through the spherical harmonic Y™. The hyper an- 
gular component ^> ang is given by 


ip* ng (1l,a) 


cosh s 




tana 


sinh 


fl(7e) (f- a ) 


(46) 


Here, we use the following hyperspherical coordinates 
1Z = \/ x 1 + y 2 = \/x 2 + y 2 
x = /? 1//2 r with /3 = 

K+ 1 

y = r 1/2 R 


x = sin ujx + cos ujy with cot ui = 


v/ 2/v + 1 


y = cos uix — sin toy 


a = arctan 


a = arctan [ — 


The function s(lZ) is determined by 


,n 


1 + s 2 


r i/2 - = ———— tanh 

a s \ 2/ 

2 cosh slo 


sinh soj 


sin 2uj cosh(s|) s sin 2 u cosh(s|)' 

The function f(TZ) is the solution asssociated with 
the lowest eigenvalue £ trimer of the hyper-radial equa¬ 
tion 


d 2 s 2 (K) + \ 


dR 2 TZ 2 

and normalised as 


- £t- 


/m = o, 


/»oo 

/ dn\fm\ 2 = l. 
Jo 


The function Ci/R) is determined by the normalisa¬ 
tion condition 

- J (sin 2a) 2 dadQ x dQyX 

a 2 

^Fad(^’ a > V) “ V , Fad( 7 ^’ V) = 

which guarantees that 

J d 3 rd 3 R \(j>(r, R)\ 2 

= ^ J lZ 5 dlZ(sin2a) 2 dadQ x dQ y R)\ 2 
= 1. 


B. Scattering of two trimers 

Trimers in the same rotational state i are identi¬ 
cal fermions and therefore scatter only in the p wave 
channel at low-energy. At sufficiently low energy, this 
p- wave scattering is negligible with respect to the s- 
wave scattering between trimers in different rotational 


























states. For this reason, we focus on the latter in 
this paper. There are three possible pairs of differ¬ 
ent rotational states, {—1,0}, {0,1}, and {1, —1}, and 
they all lead to the same scattering length, because of 
the SU(3) symmetry of this system. However, this 
symmetry is artificially broken by the single-channel 
RGM, if rotational states are given by the usual spher¬ 
ical harmonics. The different values of scattering 
lengths for the different pairs of states would thus 
give an indication of the error of the single-channel 
RGM approximation. However, a more serious issue is 
that spherical harmonics are complex-valued and the 
RGM does not ensure the scattering length to be real. 
To circumvent this problem, we consider an alterna¬ 
tive basis for rotational states, which is the xyz basis 
formed by rotational states with angular momentum 
projection zero on the three axes of space. Namely, 


Y? 


n 1 - y} 

V2 


-1 


Y? 


Y^+Yf 1 
iV 2 ’ 


Y i = Y®. 


In an exact calculation, it makes no difference whether 
one uses the usual spherical harmonics or the xyz ba¬ 
sis, but in the case of the RGM, the xyz basis ensures 
the results to be real, since the wave functions in this 
basis are all real, and restores the SU(3) symmetry 
as well. This is evident if one observes that the three 
pairs {xy}, {yz}, and {zx} can be transformed into 
each other by a rotation in space. 


To apply the RGM to this scattering problem, we 
set (f>A = 4>x and cf>B = 4> y (i.e. the two clusters are 
two trimers in rotational state x and y). There are 
twelve possible permutations of identical fermions be¬ 
tween the two trimers, as shown in Fig. [7] From this 
we obtain the expressions for the direct potential, the 
exchange potentials and the exchange kernel, which 
are given respectively by the following nine, six, and 
nine-dimensional integrals: 


with 


Vb(s) 



d 3 Rd 3 r'd 3 R 1 (\Mr-, R)M r ',&) f + 10*O', R’W y (r + M)\ 2 ) 


r± = - 


K + 1 

2k -{-1 


R± 


K + 1 

- S + 

K 


K + 1 

2k + 1 


R' - 



K 


(47) 


(48) 


I4xi(s,s')= g\jd 3 rd 3 R ( (<j> x (Ri)<f>y{ji) T <t> y (Ri)<f> x (R)J (</> x (R-2)<t> v CR-3) =F 'OC^OO'^)) (49) 

+4 {MR'MK) T MR'i)MK))* (mk'MK) t 

n 

-4 =f tvWWxmy (mk'dmkz) t) 

hi / 

V BX2 (s, s') = gX j d 3 rd 3 R ^ [^(Tl^yiRi) =F <t>y(m)^<c(Ri)) ((j>x(R-)4>y(m) T <!>y {R-)^x{R-h)T (50) 
+4 t (MK)MK' 5 ) t 0, W0*(^5))* 

hi 

-4 (0x(^ , /)0y m T cf> y (n")MK))* ^ x (Rj^y{R .") =F )■ 

hi / 


K (s, s') = 


\ I d 3 rd 3 Rd 3 r'^M^3)M^) ~F 0y(^ 3 )0*(^ 4 )) (<j> x (R)<j> v (r', R 5 ) T 0 y (^)0*(r-', R 5 ))* 

+ 4 (0x(7^)0O^) T MK)MK)) (0*W0 B (r / , R' 5 ) =F 0vW0-(r , ) H'))* ). (51) 

hi / 


In these expressions, we have set 


A = 2 


(2k + l) 6 
(2k) 3 


i (R) = ~ lim 

r— >0 or 


0 ,R) 


and R-i stands for (r,;, R, ). These variables are given 
explicitly in terms of r and R in the Appendix. The 
sign =F in Eqs. (49]|50l is — for even scattering waves, 
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Figure 7: Schematic representation of the first six permutations of identical fermions between two trimers. The next six 
permutations are obtained by performing the permutation {1 0 4, 2 5, 3 6} (which produces a minus sign for the 

corresponding terms in the wave function). 


z 


'Tl'jl'nesr—trimar nntantial \h* It? n~prW 



Figure 8: Integrated potential in the RGM1 - see Eq. {531 
- between a trimer A in rotational state x and a trimer B 
in rotational state y, as a function of the distance s and 
angle 6 (in the spherical coordinates represented in the 
inset) between the centres of mass of the two trimers. 


and + for odd scattering waves. Since we are inter¬ 
ested in s-wave scattering, only even waves are in¬ 
volved due to the conservation of parity, and thus 
=F = — in this case. Note that the asterisk in Eqs. (49 
511 still denotes the complex conjugate, although in 
our calculations all wave functions are real. To com¬ 
pute these high-dimensional integrals, we resort to 
Monte Carlo integration using importance sampling. 

The total potential (sum of direct and exchange po¬ 
tentials) 

Vi(s, s') = V d (s)S 3 (s - s') + V EX1 (s,s') (52) 

is anisotropic, due to the anisotropy of the trimers. 
To visualise this anisotropy, we plot in Fig. [8] the in¬ 
tegrated potential 

Vi integrated (^) = V D (s) + j dVl4 xl (s, s') (53) 


partial waves £ = 0 and l = 2 to get converged results. 
Fig- U also indicates that the potential is repulsive. 
This fact is confirmed by the numerical calculation of 


the potential in each partial wave given by Eqs. (23 


25). As a result, the trimer-trimer s-wave scattering 


is positive. 

We have found that the exchange potentials are to 
a good approximation local potentials. In view of 
the previous results for dimers, we substitute the ex¬ 
change potential by their local approximation given 
by Eqs. (251 and (261 and neglect the exchange ker¬ 
nel, which is costly to evaluate. We therefore use the 
local RGM1 and RGM2 equations, Eqs. ( 27p8 1. 


The resulting trimer-trimer s-wave scattering 
length is plotted in Fig. [9] as a function of the mass 
ratio k. The results are similar to the dimer-dimer 
case. As in the dimer-dimer case, the local RGM1 
and RGM2 results are very close, suggesting that the 
local approximation is enough to reproduce the RGM, 
and the contribution from the exchange of particles is 
small compared to the direct contribution. However, 
unlike the dimer-dimer case, the scattering length de¬ 
creases with the mass ratio. This is due to the fact 
that the binding energy of the trimers increases, and 
thus their size reduces, as the mass ratio increases. 



as a function of the distance s and angle 8 of the spher¬ 
ical coordinates (s, 6, p) - note that the potential does 
not depend on ip by rotational symmetry along the 2 
axis. Fig.[8]shows that the anisotropy of the potential 
is moderate. As a result, we only need to consider the 


Figure 9: Trimer-trimer scattering length att in units of 
a as a function of the mass ratio, calculated by the local 
RGM1 and RGM2. The shaded area indicates the statis¬ 
tical uncertainty due to the Monte Carlo integration used 
to calculate the potentials. 
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The decrease of the scattering length is therefore a 
consequence of the decrease of the scattering cross sec¬ 
tion, due to the decreased size of the trimers. 

V. CONCLUSION 

We have applied the resonating group method to 
the scattering of universal clusters which are described 
by zero-range interactions. We have found that the 
single-channel RGM is relevant to clusters made of 
fermions. It reproduces qualitatively, and in some 
limits quantitatively, the exact results for scattering 
involving universal dimers. We have also applied 
the single-channel RGM to the scattering of univer¬ 
sal trimers. It is found to be similar to the scattering 
of dimers: there is little contribution from the ex¬ 
change of particles and the effective interaction is re¬ 
pulsive, unlike the scattering of a fermion and a dimer, 
where exchange is dominant and produces attraction 
related to the Efimov effect. As a consequence, we ob¬ 
tain a positive trimer-trimer s-wave scattering length. 
This result has implications for the nature and stabil¬ 
ity of the ground state of a mixture of heavy and light 
fermions which are to be discussed in a separate work. 

The validity and accuracy of the present RGM cal¬ 


culations are limited by the single-channel approxi¬ 
mation. In particular, it is likely that trimers excite 
into the nearby dimer-particle continuum during their 
collision, by analogy with nuclear systems where ex¬ 
cited channels play an important role [42| . Including 
these extra channels, i.e. states of the form Eq. Q 
constructed with other eigenstates of the n-body and 
N — ?r-body subsystems, should converge to the exact 
results. It remains however numerically challenging to 
go beyond the single-channel approximation for clus¬ 
ters of more than two particles. As it stands, the 
single-channel RGM can already give useful insights 
on the interactions between universal clusters. It 
could be used to further investigate similar problems, 
such as scattering of dimers and trimers, involving un¬ 
polarised fermions or three-component fermions. 
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VI. APPENDIX 


Here we give the explicit expressions for the vari¬ 
ables appearing in Eqs. (49) and (501. 

For Eq. (|49l we have: 


„ (2k+1)(2k-1) (2k + l) 2 , 2k +1 

Ri = - 7 T - -a--—-——s H- —r + R 


2k 


2k 


k + 1 


T-2 = (k+-)(s'-s) 

„ (1 + 2k)(2k 2 - 1) (1 + 2k)(2k 2 + 2k+1) . 

-fi2 = -r—--t- S --——--- S 


2k(1 + k) 

1 T 2k 

+ -- r + R 

1 + K 

r 3 = (K+i)(s-s') + r 

1 H - 2k . .. 

R3 = 2(TT V) {s ~ s) + r 


2k(1 + k) 


. 2k -b 1 (2k + l) 2 # 2k -\- 1 

R-\ == ———~—S — -—— - S — —7--—77 -V ~b Jll 


2k 2 


2k 2 


k(k + 1) 


r 2 = (l + -)(s-s') 


r' 9 = 


(1 + k — k 2 )(1 + 2k) (1 + 2k)(k 2 + 3k 4-1) , 


2k 2 (1 + k) 

1 + 2k 

r + R 


2/t 2 (l + k) 


(1 + k)k 
1 


r 3 = (! + -s)+r 

* = 2 lITd (s - s ' , + H = H3 


„ _ ^ 1 +2k _ (1 + 2k) 2 , 

1 “ 2k(1 + k) 2k(1 + k) 

1+2k 1 + K — K 2 

(1 + At) 2 At(l + At) 

1 D 

r — R 


1 + At 


K = - 1+iT- (s + s 7 ) 


2 k 
1 + 2 At 

(1 + At) 2 

1 


1 + At — At 2 

r - 7 ^ R 


t(l + At) 


= v 1 + Yk ) ^ _ + r - r 3 

„ _ 1 + 2 At 

3 “ 2(1 + At) 


(s — s ) + R = R^ = R 3 
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Note that R" = R' 3 = R% and r" 


= r 


Additionally, for Eq. (50) we have: 


3 - 


- ~ fi ~TTT r + (1 + ^ )(s + s,) 

r 5 = -(«+ 2 X s - s ') 

R 5 = -Vl±V r _R + (2 ^ + 1)2 s /+ ( 2k + !) 

k +1 2k(k+1) 2k(k +1) 






.. 1A 2k 1 . . /\ 

R -(iT^^ (1 + s; )(8+s> 
-< 1 + ^- 8 '> 


(2ac + 1) (2/c + l) 2 , (2/c + 1) 

-7-r T* — -7-— S — 7-— S 

k(k+ 1) 2k(k + 1) 2k(k +1) 


And for Eq. (511 we have: 


t 4 = r' + (k + ^)(s - s') 
r 4 = r ’' + ( 1 +7rT)( s - s ') 












